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We use a random gap model to describe a metal-insulator transition in three-dimensional semicon- 
ductors due to doping and find a conventional phase transition, where the effective scattering rate 
is the order parameter. Spontaneous symmetry breaking results in metallic behavior, whereas the 
insulating regime is characterized by the absence of spontaneous symmetry breaking. The transition 
is continuous for the average conductivity with critical exponent equal to 1. Away from the critical 
point the exponent is roughly 0.6, which may explain experimental observations of a crossover of 
the exponent from 1 to 0.5 by going away from the critical point. 

PACS numbers: 72.20.-i, 73.40.Qv, 72.20.-i 



I. INTRODUCTION 

Undoped semiconductors have a small gap between the valence and the conduction band, typically of 
the order of 0.2. ..1.2 eV [l|. This gap is strongly affected by doping. In particular, sufficiently strong 
doping closes the gap such that a metallic phase appears. A classical example for this type of metal- 
insulator transition is doped silicon, where typical dopants are phosphorus (Si:P) or boron (Si:B) 0^. 
Disorder plays a crucial role in these materials due to the inhomogeneous distribution of the dopants. 
This suggested that Anderson localization must play a crucial role in these systems, where the quantum 
states would undergo a transition from extended to localized states for increasing disorder. This transition 
should be reflected in the transport properties, where extended states lead to a metal and localized states 
to an insulator at vanishing temperatures. 

Measurements of the conductivity <r(N) as a function of doping concentration N in Si:P at low tem- 
peratures has indeed revealed a critical behavior. Above a critical concentration N c a power law was 
found 

a(N) ~ a (N/N c - 1)" (N > N c ) 

and a vanishing conductivity for N < N c . The exponent \i was determined as fi « 0.5 for some exper- 
iments j whereas a crossover from fi fa 0.5 at some distance from the critical point to [i « 1 in a 
vicinity very close to N c was observed in other experiments @ . 

Although the picture of an Anderson transition is quite appealing, an alternative description can 
be provided by a random gap model. The idea is that the dopants create energy levels inside the 
semiconductor gap. These levels are associated with states that can overlap with the states in the 
semiconductor bands and eventually fill the semiconductor gap by forming extended states. The effect 
can be described by a random distribution of local gaps. Then the locally filled gaps can be distributed 
over the entire system and form eventually, after sufficient doping, a conducting "network". This is 
associated with a second-order phase transition which will be described in this articles. The transition is 
distinguished from the Anderson transition by the fact that the metallic phase appears at strong disorder 
(i.e. high dopant concentration) and the insulating phase at weak disorder. This does not rule out an 
Anderson transition if we increase the disorder inside the metallic regime. However, in realistic systems 
it is more likely to see the transition caused by the random gap than the more sophisticated Anderson 
transition. 



II. MODEL AND SYMMETRIES 



We consider a two-band model with a symmetric Hamiltonian. This can be expressed in terms of Pauli 
matrices Cj (j — 0, 3). A simple case is 

H = hiai + h 3 a 3 (1) 

with symmetric matrices hi, h 3 in three-dimensional (real) space. To be more specific, we can choose 
the Fourier components hi — kj \j2m with k = \Jk\ + k\ + fc| . For a uniform gap h 3 = A we obtain 
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the dispersion = ±^Jk 2 /2m + A 2 . For small wavevector k we can approximate the spectrum as 
i?k ~ ±[A + k 2 /2mA + o(/c 4 )], where the gap A renormalizes the quasiparticle mass m. Subsequently 
we will consider a random gap ft. 3 to describe the effect of an inhomogeneous distribution of dopants. 
The one-particle Hamiltonian H is invariant under an Abelian chiral transformation: 

e°" 72 He aa2 = H . (2) 

In order to reveal the relevant symmetry for transport in this system, we construct the two-body Hamil- 
tonian 



H 



H 
H 



where the upper block H acts on bosons, the lower block H on fermions. The reason for introducing 
this two-body Hamiltonian is that we can transform the distribution of the random Hamiltonian H into 
a distribution of the Green's function G(z) = (H — z)^ 1 [7, 8], which is often called a supersymmetric 
representation of the Green's function Q . 
Next we introduce the transformation matrix 

U=(l "?) (3) 



ip a 2 

and obtain the anti-commutator relation 

[H,U] + = 0. (4) 
This implies the non- Abelian chiral symmetry 

e°He (l = H , (5) 

which is an extension of the Abelian symmetry The Green's function G(z) does not obey this 

symmetry for z =/= 0. Therefore, z plays here the role of a symmetry-breaking field. An interesting limit 
is z — > 0, which we will study in the next section. 

Now we consider the case of a random gap with mean (/i3, r ) = A and variance (hs T hs T t) — A 2 = gS r r i 
and its effect on the average conductivity at frequency u. The conductivity is obtained from the Kubo 
formula as [§, LLOl 



a,, : : -— u? limiZe \ J2 r lTr2 [{G 0r {cj/2 + it)G^(-u/2 - ie))} \ , G{z) = (H - z)- 1 . (6) 



In particular, we are interested in the DC limit lo — > 0. This limit restores the chiral symmetry of H 
in ([5]) for the Green's functions. However, the symmetry can be spontaneously broken now. Since it is 
a continuous symmetry, this creates a massless mode, which represents fluctuations on arbitrarily large 
length scales. 

Here it should be noticed that <J2(H + z)~ 1 (T2 — ~(H ~ z)^ 1 . This has the consequence that the product 
in ((6]) reads Go r (z)G r o(— z) = (H — z)q t (H + z)^ 1 = — (H — z)^a2(H — z)~^a2 such that elements of 
G(u>/2 + ie) are sufficient to express the conductivity. 

III. SELF-CONSISTENT APPROXIMATION 

The self-consistent Born approximation of the average one-particle Green's function reads 

(G(*))«G (z + M7), G (z) = ((H) - z)- 1 , (7) 

where the self-energy r\ is a scattering rate, which is determined by the self-consistent equation 177 = 
Go,o( z + ill) E2- This reads in our case 



ir) = j(z + iff) 



a , ( a + X 
A - - log 



2 V a — A 



(7 = g/27r 2 , a = ^(z + ir)) 2 - A 2 /4) 
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FIG. 1: Phase diagram of the metal-insulator transition of the three-dimensional random gap model from Eq. 
([9]), where disorder is the parameter 7 and the average gap is A. 



and for z = this simplifies to 77 = r\l with 



I = 7 [A - P arctan(A/£)] , (3 = y/?f + A 2 /4 . 
In this case there arc two solutions of the self-consistent equation, namely 77 = and rj ^ with 



1 



7 



A - /3 arctan(A//3) 



(8) 



A nonzero rj reflects spontaneous symmetry breaking with respect to ([5]). Such a solution exists for ([5)1 
only at sufficiently large 7. Moreover, 77 vanishes continuously as we reduce 7. Then there is a phase 
boundary which separates the symmetric and the symmetry-broken regime: 



7(A) = 



2A - Aarctan(2A/A) 
which is plotted in Fig. [TJ The average density of states then reads 

p(E) = -5- lim Im {Tr 2 [(G rr {E + ie))}} w ^- lim Im {Tr 2 G , {E + ie + irj)} 



(9) 



Im < (E + irf) 



a , /a + A" 
A -2 l0g l^A 



) }, a= j{E + irjf- A 2 /4) 



(10) 



As a qualitative picture the average density of states is plotted for a fixed 77 in Fig. [5] 

IV. SCALING RELATION OF THE AVERAGE TWO-PARTICLE GREEN'S FUNCTION 

A common approximation for the conductivity in © is the factorization of the averaged product of 
Green's functions as USM 



E^ Tr 2[(Go r (2/)G r0 (-y))]«^r2Tr 2 [(G 0r ( 2 /))(G r o(- 2 /))] (y = w/2 + ie) , (11) 
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FIG. 2: Average density of states of the three-dimensional random gap model for fixed r\ — 0.04 and average gap 
A = 0.4 (full curve) and A = 0.8 (dashed curve). 



which can be combined with the self-consistent Born approximation in Eq. (JT]) to obtain 

E *i Tr * [(Gor(y)) (G rt ("»)>] « E r ' Tr 2 t G 0,r(y + S»?)Go,-r(-y " ^)] • (12) 
r r 

For the expression ^ this approximation leads to the Boltzmann (or Drude) conductivity, which reads 
in our specific case 

e 2 oj 2 f x A 2 1 A - z 2 
akk *2h^ J (A 2 /4-z 2 + fc 2 )3 fc dk {Z = " /2 + ■ (13) 

Thus the conductivity vanishes in the DC limit to — > for r\ > 0. The reason is that the self-consistent 
Born approximation creates the Green's function Go, r (?/ + if]), which decays exponentially on the scale 
1/rj. Consequently, the sum over the real space on the right-hand side of Eq. (fT2")l is finite. 

A more careful inspection indicates that the averaged product of Green's function on the left-hand side 
of Eq. (fTTj) decays according to a power law as a consequence of the massless fluctuations around the 
spontaneous symmetry breaking solution 77 ^ Q • We can perform the integration with respect to these 
fluctuations and obtain the relation [8( 

E r l Tr2 [(Gor(y)G r0 (-y))} = f(r}/y) E r 2 k Tr 2 [G , r (y + irj)G ^ r {-y - irj)] , (14) 

r r 

where the coefficient depends on the ratio of the order parameter of spontaneous symmetry breaking 77 
and the symmetry-breaking field y: 

f{i 1 /y) = {l + n 1 /y) 2 . (15) 

This coefficients represents correlations of the Green's function fluctuations, which are negligible only for 
f(r)/y) ~ 1. In the absence of symmetry breaking is rj — and /(0) = 1. This justifies the approximation 
by Eq. (fT2"]) in the insulating regime. Moreover, f(i]/y) diverges in the presence of spontaneous symmetry 
breaking for a vanishing symmetry-breaking field y. In particular, it is finite for u> > which reflects the 
fact that fluctuations are cut-off on length scales L w ~ vp/u (vp is the Fermi velocity). 
With the scaling relation (fl4l the conductivity in Eq. © becomes in the DC limit ui — s- 

e 2 j] 2 

°kk = — -—==== . (16) 
47m ^/A 2 /4 + rj 2 
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FIG. 3: Conductivity as a function of disorder for an average gap A = 0.004 (dotted curve), A = 0.04 (full 
curve) and A = 0.4 (dashed curve). There is a metal-insulator transition at 7 » 1, at 7 » 1.03 and at 7 « 1.37, 
respectively. The second plot demonstrates the fit (dashed curve) of the conductivity (full curve) away from the 
critical point 7 C for A = 0.01 by 0.47(7 - 7c) ' 6 - 



The solution 77 of the self-consistent Eq. © is inserted into akk and the conductivity is plotted as a 
function of disorder strength 7 in Fig. [3J The conductivity vanishes linearly with decreasing disorder 
strength (i.e. with decreasing doping concentration). 

V. DISCUSSION AND CONCLUSIONS 

Our result for the DC conductivity in Eq. (|16|) . together with the solution of the order parameter 
r\ in Eq. ([8]), provides a simple description of a metal- insulator transition in doped three-dimensional 
semiconductors. The metal-insulator transition is characterized by the scattering rate 77 that vanishes 
in the insulating regime. Such a behavior is not an Anderson transition, since the latter would have 
a scattering rate i) ^ on both sides of the transition [l2j ■ Even more important is the change of the 
coefficient f(rj/y): It is always 1 in the insulating regime and infinite in the metallic regime. This quantity 
describes the correlations of the Green's function fluctuations in the relation (ITU) . 



G 



There is a linear behavior near the metal-insulator transition and a crossover to a non-critical power 
law, as depicted in Fig. [3J For the linear part the slope of the conductivity is quite robust with respect 
to the average gap A (cf. first plot in Fig. [3]). Away from the transition point a negative curvature 
appears though, which can be fitted by a power law with exponent /i s» 0.6 (cf. second plot in Fig. [3j. 
The change of exponents can be related to the discussion in Refs. 0, [l3| about a crossover of exponents 
in Si:P from /j, ss 1 very close to the critical point N c to /i ~ 0.5 further away from N c . Rosenbaum et 
al. have found that the conductivity close to the critical point varies from sample to sample 5]. This 
indicates strong conductivity fluctuations, which may also exist in our random gap model, as indicated 
by the strong fluctuations of the Green's functions due to the large values of f(rj/y). 

The conductivity always increases with -q. This may be an artifact of our approximation because 
we have neglected (massive) symmetry-breaking modes in the derivation of the scaling relation Eq. 
(fl"4|) . Symmetry-breaking modes can become unstable for sufficiently strong disorder 7 > 7, > -f c . 
These unstable modes terminate the validity of our approach, as we have seen for two-dimensional Dirac 
fermions with random gap [7j. Such an instability might be linked to Anderson localization. If this is 
true, it would be of interest to increase disorder in the metallic phase of doped semiconductors to observe 
a genuine Anderson transition. Another interesting aspect in these materials is the insulating phase, 
where magnetic correlations play a crucial role [Til fl5| . This would require an extension of the random 
gap model by including spin degrees of freedom. 
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